Phase diagram of the two-dimensional t— J model at low doping 
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The phase diagram of the planar t-J model at small hole doping is investigated by finite size scal- 
ing of exact diagonalisation data of \/~N x y/~N clusters (N < 26). Hole-droplet binding energies, 
compressibility and static spin and charge correlations are calculated. Short range antiferromag- 
netic correlations can produce attractive forces between holes leading to a very rich phase diagram 
including a liquid of d-wave hole pairs (for J ft > 0.2), a liquid of hole droplets (quartets) for larger 
J/t ratios (J/t > 0.5) and, at even larger coupling J/t, an instability towards phase separation. 

PACS numbers: 74.72.-h, 71.27.+a, 71.55.-i 



Studying the behavior of holes in two-dimcnsionnal 
(2D) antiferromagnets is crutial to understand the ori- 
gin of pairing in the high-T c cuprate superconductors. In 
these materials chemical substitutions in the parent sto- 
chiometric compound lead to injecting mobile holes in the 
CuOi antiferromagnetic planes. Besides transport prop- 
erties these holes will also drastically affect the antifer- 
romagnetic (AF) correlations in the planes. On the the- 
oretical side, motion of holes in antiferromagitets can be 
simply described by the so-called t-J modefiEl, a strong 
coupling version of the .well-known Hubbard model. Pre- 
vious numerical studiesu have given reliable informations, 
specially in the limit of a single hole. It is believed that 
holes behave like quasiparticles, at least at the bottom 
of the coherent band, although spin fluctuations strongly 
enhance theix effective masses and reduce their quasipar- 
ticle weightsa 

Finite size scaling analysis becomes easier at commen- 
surate densities such as n = 1/2 (quarter filling). In 
this case, exact diagonalisations (ED) studies of the t- 
J modelu have suggested the existence of superconduct- 
ing correlations in the vicinity of the phase separation 
phaseu. This regime is however quite far from the exper- 
imental situation. 

At small but finite doping (eg electron density n ~ 0.8- 
0.9) fewer theoretical results are known. However it is b&, 
lieved that this class of models reproduces successfullyEI 
the large Fermi surface observed in angular resolved pho- 
toemission studies in the metallic phase of doped high-T c 
materials. In addition, possible observation of shadow 
bands due to strong short range antiferromagnetic cor- 
relations has been, suggested in both experimentafl or 
theoretical studies^. 

The magnetic coupling J can generate an effective cou- 
pling between holes. This is particularly clear in the 
(unphysical) large J/t regime where the magnetic en- 
ergy cost is minimized by having holes sitting on near- 
est neighbor sites. In this regime, the uniform state be- 
comes in fact unstable towards_a phase separated stateQ. 
High temperature expansionsL9 also predict phase sepa- 



ration for J/t > 1. However, small cluster calculations 
have shown that, for smaller and mcpejiealistic J/t ra- 
tios individual pairs could be stablalil't3. Preliminary 
resultal3 state that larger clusters of holes could also form 
in the intermediate parameter range. Other possible can- 
ditates,i|ii,this parameter regime are non-uniform striped 
phasetM 

In this letter, more insights into the nature of the phase 
diagram at small doping are obtained from a detailed 
numerical study. Indeed, since analytic perturbation 
treatments are poorly controled in the relevant physical 
regime, exact diagonalizatkms of small 2D square clusters 
by the Lanczos algorithmic were performed!!!!. Studies 
in the regime of small finite hole densities are delicate 
since only different discrete values of the densities can be 
achieved on different clusters and interpolations between 
them become then necessary. First, finite size scaling of 
binding energies of n-hole clusters provides indication on 
the stability of liquids of pairs or droplets in the van- 
ishing hole concentration regime. Hole-hole correlations 
obtained for n ~ 0.85 also confirms the stability of pairs 
at finite density even at small J/t ratios. The compress- 
ibility for arbitrary hole densities (n < 0.8) and various 
system sizes is calculated to perform an extrapolation 
to the thermodynamic limit. The domain of the phase 
separated region is then estimated. We also discuss the 
behavior of the static spin and charge structure factors 
at intermediate density n ~ 0.85. 

The t-J model defined on a square lattice reads, 



H = t £ 4 

x,y N.N. 
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where c£ a and S x are creation and spin operators at site 
x. Ground state energies (GS) and equal-time correlation 
functions in the GS are obtained on small \J~N x \f~N N- 
site clusters at low hole densities by the Lanczos method. 
Typically N=18, 20 and 26. 

Let us first consider a fixed finite number of holes 
Nh = n (Nh = 2, 4) on various clusters of increasing 
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sizes. These holes will form a n-particle bound state if 
the binding energy A„ = E h , n + E h , - 2E h ^ n / 2 , where 
Eh,n is the GS total energy for a system with Nh = n, 
converges towards a negative value in the limit of infinite 
system size. Strictly speaking such quantities give indi- 
cations about the stability of n-particle boundstates only 
in the limit of vanishing hole density. 

A simple broken-bond counting argument shows that, 
at large J/t, two holes injected in the AF can minimize 
the local magnetic energy by, forming a bound state. Fi- 
nite size scaling have shownEJ that this picture is actually 
correct even down to J = Jb.2 — 0.2. A 2 calculated for 
N=26 is shown in Fig. [I] and becomes negative when the 
paired state is stable. The hole- hole pair has a d x 2_ y 2 
orbital symmetry. Fig. ||(a) shows the hole-hole density 
correlations of the pair Cu{r) = Nh ^ h _^ (n/ I (r)n/ t (0)) 
for all possible distances (r ^ 0) between the holes com- 
patible with the cluster shape. Note that the normal- 
ization factor is chosen so that X) r ^o^' 1 ( r ) = 1- Cor- 
relations at the intermediate distance of \[2 ie when the 
holes stay across the diagonal of a plaquette on the same 
sublattice are dominanttJ. This singlet hole pair can in 
fact be viewed as a combination of a triplet pair with 
a nearby spin triplet excitationE3. The two-hole X1S ex- 
hibits flux quantization in a ring in units of /ic/2al3, also 
a signature of a paired state. Lastly, pair formation is 
also consistent with the observation that the dynamical 
response of the d-wave pair creation operator exhibits on 
small clusters a sharn_<5-peak of weight Z 2 h oc J/t down 
to small values of J/tO. 

We now consider the possibility of larger droplets of 
holes. Indeed, above Jb,2 residual interactions between 
pairs might be sufficiently attractive to stabilize e.g. hole 
quartets. GS energies of 4-holes can be calculated for 
three possible orbital symmetries (s, p or d-wave) of the 
wave function (with zero total momentum). The lowest 
energy is obtained in the s-wave channel. In order to 
estimate the onset of clustering A4 has been evaluated 
on various lattices of size up to N=26 and the data are 
displayed in Fig. [l| Note that the critical value J B ,4, at 
which A4 changes sign depends weakly on the system 
size while the slope \dA4/dJ\ at this point decreases for 
increasing size. 

Fig. [I] then strongly suggests that when J/t S 
[J/t\ B ,2,J/t\ B ,i] with J/t\ B .2 ~ 0.2 and J/t\ BA ~ 0.5 
individual hole pairs exist without forming larger clus- 
ters. This results based on a scaling of A„ is, strictly 
speaking, only valid in the limit of vanishing hole den- 
sity. In order to investigate the stability of the pairs at 
finite doping we have calculated the hole-hole correla- 
tions Ch(r) in a 26-site cluster at density n ~ 0.85 (i.e. 
with 4 holes). The results shown in Fig. ||(b) for various 
separations |r| reveal dominant correlations at distance 
y/2 as in the case of a single pair (see Fig. |](a)). Note 
also that the density correlations at the largest distances 
available in the cluster remain always significant in this 
parameter regime which is consistent with the existence 



of separate pairs in the cluster. 

Formation of droplets should not be confused with 
phase separation (PS) between hole-free and hole-rich 
phases even though both phenomena have the same mi- 
croscopic origin. The issue of phase separation can be 
addressed by studying the inverse compressibility kT 1 de- 
fined by k _1 = d 2 (Eh, n /N) / dn\. In a uniform system 
k _1 is finite and positive. On finite clusters kT x < sig- 
nals the instability of the homogeneous phase. However, 
we stress that a rigourous determination of the PS region 
can only be achieved by a finite size scaling at constant 
hole density. Such a study is attempted in Figs. ||(a)(b) 
showing the GS energy per site vs hole density nh for 
two different system sizes. The thermodynamic limit is 
obtained in two steps; (i) an interpolation between data 
points at constant cluster size and (ii) an extrapolation at 
constant hole density assuming N- 5 ' 2 finite size correc- 
tions. We have checked by considering other system sizes 
(data not shown for clarity) that such a TV -3 / 2 behavior 
is actually very well satisfied provided that J > n^t. For 
J > 1 it is clear that the PS region extends from nh = 
to nh — 0.12. For smaller values of J becomes 
quickly very small and ws-can assume that there is no 
sign of PS for J/t < 0.75Ej consistently with the results 
obtained by high-temperature expansions!!!! Our results 
then suggest that a liquid of hole quartets are stable in a 
small region of the phase diagram, as a precursor of the 
PS instability line. This should be contrasted to the small 
electron density, case where the gas of electron quartets 
is never stablecJ. 

We finish this study by the investigation of the 
static spin and charge structure factors defined 
by 5(q) = |L(^W5 z (0))^ r and N(q) = 
jfJ2 T ((nh(r)nk.(0)) — n 2 i )e l ' i ' r respectively. The data ob- 
tained for N=26 with a density of n ~ 0.85 (four holes) 
are shown in Figs. |4](a) and (b). A smooth interpolation 
between the discrete q-points of the reciprocal lattice of 
the 26-site cluster has been performed assuming that the 
correlations in real space remain small at distances larger 
than the cluster size. S(q) in Fig. |](a) shows a pro- 
nounced peak at (tt, w) even for small values of J. This 
indicates that large anti-ferromagnetic spin correlations 
((,af ~ 3) still survive for hole doping as large as 15%. 
iv(q) shown in Fig. [|(b) exhibits along the T-M line a 
behavior very similar to non-interacting spinless fermions 
with nearest-neighbor hopping (dotted line). However, a 
clear dip is observed at X. This behavior cannot be ex- 
plained by a simple Fermi surface effect. For example, 
a different tight-binding spinless model whose dispersion 
has a minimum at £ (momentum (w/2,tt/2)) would give 
much more structure than observed (other dotted line). 
We interpret the dip at X as the signature of strong 
short-range correlations between holes characteristic of 
the paired state. 

We conclude this paper by suggesting a possible phase 
diagram in Fig. || based on the results discussed above. 
When J exceeds some critical values Jb,2 and Jb,4 holes 
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injected into the antiferromagnetic phase form a liquid of 
d-wave pairs and a liquid of quartets (D) respectively. At 
larger J/t ratios the t-J model phase separates (PS). Also 
note the existence of a ferromagnetic region (F) at vecj; 
small J-.&S predicted by high temperature expansionsEHI 
or EDE3. Two crutial issues still remain to be addressed 
namely the exact nature of the normal paramagnetic 
phase (P) and possible pair-pair correlations (supercon- 
ductivity) in the pair liquid phase. 
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FIG. 1. Binding energy A4 vs J for clusters of 16, 20 and 
26 sites. The 26-site cluster hole-hole binding energy A2 is 
also indicated by open stars. 

FIG. 2. Hole-hole correlation function for various hole sep- 
aration vs J/t obtained on a 26-site cluster with two (a) and 
four (b) holes. The various symbols associated to the allowed 
distances are indicated on the figure. 

FIG. 3. (a) energy per site vs hole density for 18 and 26 
site clusters and various J/t values (indicated on the plot). 
Continuous lines are interpolations between data points cor- 
responding to the same system sizes. Extrapolation to the 
thermodynamic limit are indicated by dotted lines, (b) data 
for J/t=l only shown on an enlarged scale. 

FIG. 4. Static spin (a) and charge (b) structure factors 
S(q) and iV(q) along symmetry lines of the Brillouin zone 
for various J/t ratios as indicated on the plot. V, M and X 
corresponds to (0,0), (7r,7r) and (71", 0) respectively. In (b) 
the dotted lines are obtained assuming the Nh holes behave 
as non-interacting spinless fermions with various dispersion 
relations (see text). 



FIG. 5. Schematic picture of the J/t - n phase diagram. 
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